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Abstract 

The effect of an ac modulation on a 2-level or 3-level system is studied theoretically. The 
absorption spectrum is calculated by solving the optical Bloch equations and is interpreted by 
reference to the Floquet quasienergy spectrum. The dependence of the absorption sidebands on 
the intensity of the coupling laser field in Rydberg systems submitted to a radio-frequency (rf ) field 
is analysed in detail. It is shown that for sufficiently strong coupling fields additional sidebands 
appear in the probe absorption spectrum in ladder 3-level systems. These additional sidebands are 
induced by the coupling of the intermediate state to the Floquet manifold spawned by the upper 
state under rf modulation. 
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I. INTRODUCTION 



There is currently much interest for electromagnetically induced transparency (EIT) in 
Rydberg atoms g in particular for the cooperative atom-light and photon-photon dynamics 
found in Rydberg systems |2H5j and for their possible application, e.g., in quantum informa- 
tion processing in the production of correlated photon pairs [7|, and in metrology js, 9|. 
The present work is motivated by a recent experimental study of EIT in a 3-level Rydberg 
system with radio frequency (rf) modulation 10[. In this experiment, the absorption on 
the 551/2-5 P3/2 transition of 87 Rb was probed by a weak laser field in the presence of both 
an rf field and a strong laser field resonant with the transition between the intermediate 
5Ps/2(F' = 3) state and the highly polarizable ?>2Si/2{F' = 2) state. The rf field was too low 
in frequency for driving multiphoton transitions between the lower, intermediate or upper 
states and other bound states. However, it was strong enough to Stark-shift the upper state 
periodically by tens of MHz, which resulted into the EIT features in the probe absorption 
spectrum shifting and splitting into multiple EIT sidebands. 

The existence of absorption sidebands is a well known feature of ac-modulated 2-level 
systems 



111 ], including 2-level atomic or molecular systems 



12 



19| . The effect of an oscillat- 



ing electric field on molecular absorption lines with quadratic Stark shift was described by 



Townes and Merritt in 1947 



13j: Lines whose frequency width is much larger than the fre- 



quency of modulation adiabatically follow the oscillations of the field — i.e., a state of dipole 
polarizability a in a non- resonant ac field £(£) = £0 sin(o;o£) is simply displaced in energy by 
a[8,(t)} 2 /2 at time t. However, lines whose frequency width is of the same order or smaller 
than the modulation frequency split into manifolds of components displaced in energy from 
their zero-field limit by a£g/4 + 2nhuo, n = 0, ±1, ±2, . . ., and the intensity of these com- 
ponents varies with n, ujq and £ like J 2 (aHo/8huo). Analogous shifts and sidebands are 
also found in systems governed by an Hamiltonian unitarily equivalent to that governing 
ac-modulated 2-level atomic or molecular systems, such as ac-driven superconducting flux 
qubits [20|427|. 

Interactions with other fields will in general introduce additional n-dependent Stark shifts. 
For this reason, the separation between adjacent sidebands of an ac modulated 2-level sys- 
tem is not exactly 2huo and the intensities of the absorption lines do not exactly follow 
a J 2 (aEo/8hcjQ) law. The correct separations and scalings can be obtained by combining 
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Floquet theory and Van Vleck perturbation theory [25l l27|. However, these deviations from 
the results of Townes and Merritt are significant only when the interaction coupling the two 
states of the system is sufficiently strong. Perhaps for this reason, they do not appear to 



have been previously considered in an atomic and molecular physics context 

In many ac- modulated 2- level atomic systems of practical interest, one of the two states 
is much more polarizable than the other and the frequency widths of both the state and the 
coupling laser field are small compared to the modulation frequency. In these conditions, it 
may be tempting to ignore the manifold structure of the less polarizable state, retain that 
of the more polarizable state, and model the ac-modulated two-level system as a many-level 
system in which the less polarizable state interacts with a manifold of independent sideband 
states. We refer to this approximate description as the iV-level approximation. It leads to the 



results described by Townes and Merritt 



13J in the limit of a weak interaction between the 



states and makes it possible to calculate the absorption spectrum rapidly at a minimal cost. 
However, as is shown below, this approximation breaks down for strong coupling between 
the two states: In order to obtain the correct results, the sideband structure of both states 
must be properly taken into account. This can be done, e.g., by carrying out a full Floquet 



analysis of the ac-modulated 2-level system 



25 



27|, 



29 



-|32|. or by solving the optical Bloch 



equations with time-dependent modulation. The sidebands of the less polarizable state here 
originate from the coupling of this state to the Floquet manifold formed by the other state, 
not from a direct effect of the ac-modulation. Nonetheless, as is illustrated below, they may 
manifest by well defined, additional resonance structures in absorption spectra. We refer to 
these as induced sidebands. 

In this paper, we investigate the changes in the sideband structure of the absorption 
spectrum arising from the interaction with the coupling laser field, both in two-level and 
three-level systems. The absorption spectrum is calculated by solving the optical Bloch 
equations and is interpreted by reference to the Floquet quasienergy spectrum. The iV-level 
approximation is compared to the exact results throughout the paper. For simplicity, we 
first focus on ac-modulated two- level systems (Sec. HI])- The coupling with a lower state is 
then introduced. The probe absorption in an rf-modulated 3-state ladder system and the 
appearance of induced sidebands is studied in Sec. IIHI In particular, we give examples of 
sidebands shift in EIT spectra for both cold and thermal ensembles. We conclude with a 
summary of our main findings (Sec. II V]) . 
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We use atomic units throughout the rest of this paper, except where indicated otherwise. 



II. RF-DRESSED TWO-LEVEL SYSTEMS 

A. RF-dressing in the adiabatic approximation 

We first concentrate on the rf dynamics of the subsystem formed by the intermediate 
and upper states of our 3- level ladder Rydberg system, states |6) and |c). We thus assume, 
for the time being, that neither \b) nor |c) interact with the lower state of the system, \a). 
Making the rotating wave approximation for the laser field and passing to slowly-varying 
variables would then reduce the time- dependent Schrodinger equation to 



d 


'0 




dt ' 




H 




(1) 

(which we 



-iVt c /2 
-itt c /2 iA c 

if the rf field was absent. Q c and A c are respectively the Rabi frequency 
assume real) and the detuning of the control laser field coupling the states \b) and |c) in the 
ladder system: Denoting the zero-field energies of these two states by Wb and w c and the 
angular frequency of the control laser by u c , A c = u c — u cb with uj c b = w c — io&. The state 
vector of the system, is a linear superposition of the bare states \b) and |c) and the 

functions q,(t) and c c (t) are the respective probability amplitudes: 

m)) = c b (t)\b)+c c (t)\c). (2) 

We write the electric field component of the rf field as 

£(i) = £ r f sin(a; r ft). (3) 

We assume that the angular frequency u r f is so much smaller than Ub c , and so much smaller 

than the transition frequencies between \b) and |c) and any other unperturbed state, that 

| b) and |c) evolve adiabatically under the effect of this field. Hence, we repi 

dependent bare states exp(— iwbt)\b) and exp(— \w c t)\c) by the adiabatic states 

-t 



ace the time- 

nn 



\B(t)) = exp (-i j l E b [Z(t')} dt'^j |&[£(t)]> (4) 



and 



\C(t)) = exp f-i ^ E c [E(t')] dA |c[£(t)]>. (5) 



In these expressions, \b[E]) and |c[£]) denote the state vectors which develop adiabatically 
from the unperturbed states \b) and |c) when a static electric field is turned on fromn to 
£, and E b [E] and E C [E] are the corresponding eigenenergies of the Stark Hamiltonian. The 
rf field is assumed to be weak enough that, in sufficiently good approximation, E b [E(t)] = 
w b -a b [E(t)] 2 /2 and E c [E(t)\ = w c -a c [E(t)] 2 /2, and moreover that \b[E(t)]} = \b) + E(t)\b^) 
and |c[£(i)]) = \c) + E(t)\c^). The two state vectors \b^) and \c^) are defined by first order 
perturbation theory, and a b and a c are the static dipole polarizability of the states \b) and 
|c). For the applications we have in mind, \a c \ ^> \a b \, so that one can take E b [E(t)} ps w b 
and |6[£(£)]) ~ \b), and assume that only the upper state |c) is significantly perturbed by 
the rf field. Moreover, we are also considering systems for which \b) and \c^) have the same 
parity and therefore are not directly coupled to each other by the control laser field. In these 
conditions, the time-dependent Schrodinger equation for the rf-dressed system becomes 

-z'fi c /2 
-ifl c /2 i[A c — 2E C sin 2 o; r ft] 

where 

S c = a c E 2 jA. (7) 

Eq. (jSJ) is a system of linear differential equations with periodic coefficients, which is 
amenable to the Floquet description studied in the following section. The effect of relaxation, 
which is not taken into account in Eq. ([6]), will be considered in Sec. Ill CI 

B. Floquet formalism 

1. Floquet states 



d 


(J 




dt ' 




H 




(6) 



By virtue of the Floquet theorem 25|, I29l-l3l| , any solution of Eq. (EJ) can be written as 



a superposition of fundamental solutions of the form 

oo 

c b {e-t) = e~ iet cj/" 2 ™"' (8a) 

■«=— oo 
oo 

c c (e;t) = e-* ^ c c , n e~ 2m ^\ (8b) 

n=— oo 

where the quasienergy e and the coefficients c b)U and c C)U are constants. (The fundamental 
angular frequency is 2cj r f, not w r f, because the time-dependence entirely arises from the 



periodic Stark shift of the upper level, which varies at twice the rf frequency.) Replacing 
Cfi(t) and c c (t) by c&(e; t) and c c (e; t) turns Eq. (jfJJ) into a time-independent system of algebraic 
equations for the coefficients c^ n and c c n , namely 

Y Cc <" - 2Tto ^" = eCb - (9a) 
s . 

— ^(c c , n -i + c c , n+ i) = ec c , n . (9b) 

The quasienergies are those values of e for which this system has a non-trivial solution. In 
matrix form, Eqs. (l9~a|) and ( )9b1) read 

H 2F c = ec, (10) 
where c is the column vector (. . . , Cb-i, c^o, q, i+ i, . . . , c Cj _i, c C) o, c Cj+ i, . . .) T and 
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2F 



V 



2kJ r f 
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fi c /2 














•• 
















— 2u; r f • ■ 
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■ S c - A c + 2u r{ 


-S c /2 
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The possible values of the quasienergy are the eigenvalues of H2F, which we denote by e^. 
Each one corresponds to a particular dressed state of the system, i.e., to a solution of the 
Schrodinger equation of the form 



-ie k t 



E 



-2inoj r ft 



(12) 



where the state vectors |"0fc,n) are time-independent. It follows from the above that = 
Cft(efc; i)|6) +c c (efc; i)|c). We thus have two equivalent descriptions of the dressed states of the 
system, namely one in terms of the time-dependent coefficients Cb(ek, t) and c c (ek, t) and one 
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in terms of the time-independent vector l^jt)) = (\ipk,o), l"0fc,±i)> ^,±2), • • ■) T formed by the 
harmonic components of (Here and in the following, we use the notation | . . .)) to 

represent a column vector of state vectors.) The vector solution of Eq. fllOp with e = e^, 
is the representation of \^k)) i n the basis of the bare Floquet states 

\b,n)) = (\b)5 m , t = 0,Tl,T2,...) T (13) 

and 

\c,n)) = (\c)5 m ,i = 0,Tl,T2,...) T - (14) 



Thus 



= Yl (cb,-n\b,n)) + c Cin \c,-n))). (15) 



The Hamiltonian matrix H 2 f has a four-block structure. Its upper and lower diagonal 
blocks, Hb and H c respectively, are square matrices representing the Hamiltonians of the 
Floquet manifolds spawned by the bare states \b) and |c). The non-zero elements of the 
off-diagonal blocks arise from the coupling of the states \b) and |c) by the control field, while 
the off-diagonal elements — X c /2 of H c arise from the coupling of the state |c) with itself via 
the absorption or stimulated emission of two rf photons. The rf field thus mixes the different 
vectors |c, n)) with each other. By contrast, the matrix is diagonal since we neglect the 
dressing of the state \b) by the rf field. 

It is useful to diagonalize the matrix H c , which can be done by transforming the basis from 
{\b, n)), \c, n)), n = 0, ±1, ±2, . . .} to n», \d, n)), n = 0, ±1, ±2, . . .}, with \b', n)) = 
\b,n)) and 

\c,m)), (16) 



m — rv^ V / 



m=— 00 



where J p (x) is the p-th order Bessel function of the first kind 25|. Under this change of 
basis, the matrix representing the Hamiltonian becomes 
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H' 



2F 



2w rf 



-2u ri 

Q c J /2 fi c J_i/2 fi c J_ 2 /2 

Vt c J_s r \/1 f2 c< /o/2 f2 c< /_i/2 

^c^+2/2 Vt c J + i/2 Q c J /2 



^c<A)/2 Q c J + i/2 Q c J + 2/2 
£l c J—\/2 Vt c J§/2 Vt c J-\-\/2 

n c j- 2 /2 fi c J-i/2 n c Jo/2 



A c + 2cu Ti 
S c - A c 








A r - 2u 



rf 



V 



; ; ; ; ; ; ; ■.) 

(17) 

where J p , p = 0,±1,..., denotes J p (S c /2o; r f). Neither the quasienergies nor the diagonal 
elements of the Hamiltonian matrix are affected by this change of basis. We see that in 
the absence of the control field, i.e., for Q c = 0, the vectors \b',n)) and |c',m)) describe the 
Floquet sidebands of the states \b) and |c), respectively. When Q c ^ 0, the sidebands of b are 
coupled to those of c by the off-diagonal blocks of H' 2F , the strength of the coupling between 
;he n-th side band of b and the m-th side band of c being proportional to | J m _„(E c /2a; r f)| 
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2. Quasienergy spectrum 

H2F and H^p reduce to the same diagonal matrix in the absence of control field, i.e., for 
Q c = 0. It follows from Eq. (TlTl) that their spectrum is a double comb of quasienergies: 
When f2 c = 0, any quasienergy is equal either to e[ ^ or to e^l for some value of n, where 

eg = 2nu rh (18a) 
eg> = S c - A c + 2nu ri . (18b) 

The corresponding dressed states, \^k)), are either \b',n)) or |c',n)), respectively. Each of 
the quasienergies is thus associated with an energy sideband of either state \b) or state |c). 
The sidebands of \b) are located at the energies Ef^ = Wb + eg = Wb + 2nu r f, n = 0, ±1, . . ., 
and those of |c) at the energies Ej?n = Wb + w c + e c°« = w c + S c + 2noj rf , n = 0, ±1, . . . - 
see Fig. [T](a). As should be expected when Q c = 0, the sideband energies and Ei°n 
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FIG. 1: (a) The manifold structure of the dressed 2-level system. The energy separation between 
adjacent sidebands is 2w r f . The energy of state |c) is Stark shifted by the rf field by an amount S c . 
The effective detuning from resonance of the laser field coupling the two states is thus A c — S c . (b) 
The manifold structure in the iV-level approximation: the sidebands of state \b) are neglected and 
the interaction between the two manifolds reduces to the interaction of the single state \b) with 
the N states forming the c manifold. 



do not depend on the frequency of the control laser. However, the quasienergies of the 
ec°2-manifold of quasienergies change with the detuning A c , with the consequence that this 
manifold crosses the (w c -independent) e[ ^-manifold at A c = £ c + 2mu r [, m = 0, ±1, ±2, . . . 

The quasienergy spectrum keeps that double comb structure for Q c ^ 0: any quasienergy 
6^ is equal either to e& jn or to e Ci „ for some value of n, with 

eb, n = e bfi + 2nu rl , (19a) 

e c , n = e c ,o + 2nu; rf , (19b) 

e h o and e Cj0 being such that e^o — > e i°d an d e c ,o — > 4°o m the limit O c — ^ 0. Thus e^^n ~ efl and 
£c,n ~ Cc°n if the control field is sufficiently weak. However, when Q c ^ 0, the eb jn -manifold 
and the e Cjn -manifold are no longer degenerate at A c = S c + 2mu r f, m = 0, ±1, ±2, . . . 
Instead, as illustrated by Fig. [2(a), these two manifolds anticross each other in sequences of 
avoided crossings. 

The gaps and positions of these avoided crossing are determined by the off-diagonal 
blocks of the Hamiltonian matrix H' 2F . If E c = 0, i.e., if the rf field is vanishingly weak, 

9 




A c /2tt (MHz) 

FIG. 2: (color online) The quasienergy spectrum, (a) and (b), and the average population in the 
upper state at large times (t 3> 1/T C , as defined in Section Til Cp . (c), in an rf-dressed 2-level system. 
O c /2vr = 10 MHz, S c /2vr = 10 MHz, w rf /2vr = 10 MHz, r c /2vr = 1 MHz and 7 C = 0. (a) and dark 
pink curve in (c): exact results, (b) and light purple curve in (c): results obtained in the iV-level 
approximation. Maxima in the population transferred to the upper state occur at the detunings 
at which the quasienergy curves cross or anticross each other. Differences in the positions and 
amplitudes of these maxima between the exact results and the iV-level approximation are visible 
in (c). The crosses represent one of the quasienergy curves predicted by Eqs. (f25l) and ([26]) . 



J p (S c /2(jj r f) = 5 p o, and these off-diagonal blocks are the unit matrix multiplied by Q c /2. In 
this case, the system is equivalent to an infinite set of non-interacting two-level systems. 
One then has 

e b , n = ~ ± ^V A c + W c + 2nuj rU (20a) 

e c , n = ~ T ^y/Al + Ql + 2nu rU (20b) 
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where the upper signs apply for A c > and the lower signs for A c < 0. The states 
interchange character at A c = 0. However, since in the absence of the rf field Floquet 
harmonic components of different n values are not coupled to each other, the crossings 
between the b and the c manifolds at A c = 2muj T f (m ^ 0) are true crossings rather than 
avoided crossings. 

The off-diagonal elements of the off-diagonal blocks are normally non-zero if the rf field is 
not vanishingly weak, but they rapidly decrease away from the diagonal when |E c /2aj r f| <C 1 
since J p (E c /2u; r f) then goes rapidly to zero when \p\ increases. If one neglects these off- 
diagonal elements, Eqs. fl20a]) and (120bj) then generalize to e& jn ~ e£„ and e c n « e^l, with 



e (i) _ I ( € <f>) , e (o) \ ± 1 // (o) _ (o)\ 2 Q2J2 

= i (E c - A c ) ± ^(E c -A c ) 2 + n?j2 + 2nu T{ (21) 



and 



e c,n — 2 + C c,n y + 2 V \ / c 



= ^ (S c — A c ) =f ^ ^/ (E c — A c ) 2 + Q^Jq + 2no; rf . (22) 

The upper signs apply for A c > E c and the lower signs for A c < E c . The control laser field 
thus lifts the degeneracy between the states \b',n) and \c',n) at A c = E c and creates an 
energy gap approximately equal to fl c J (E c /2c<; r f) between e^ n and e C)U . In Fig. EJja), this 
interaction results in conspicuous avoided crossings at A c /(27r) ps 10 MHz. 

The two manifolds formed by the quasienergy curves e^ n cross each other at the detunings 
at which efl = e^ ±m , m = 1, 2, . . ., i.e., at A c = A Ci±m with 



A Ci±m = E c ± ^4m^ { - ttPl (23) 

Therefore these true crossings do not occur exactly at A c = E c ± 2mu r f, as would be the 
case if Q c = 0, but instead at detunings shifted towards the avoided crossings at A c = E c . 
Expanding the square root function to lowest order in Q c , Eq. (J23l) gives 



A c , ±m « E c ± 2mw ri =f ^ c 2 J 2 /(4mu rf ). (24) 
This result can also be obtained using perturbation theory. 
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The degeneracy of the eft and eft ±m quasienergies at A Cj±m is lifted by the off-diagonal 
terms of the off-diagonal blocks of the Hamiltonian matrix H' 2F . For instance, the avoided 
crossings visible at A c /(27r) « —10 MHz and 30 MHz in Fig. [2(a) originate from the terms 
O c J±i in these off-diagonal blocks, which couple the \U,n)) states to the |c',±l)) states. 
Treating each pair of intersecting quasienergies as an isolated two-level system in the vicinity 
of their intersection and neglecting the terms in Q c J± m with m > 1 in H 2F yields the improved 
approximation e M » eft and e c ,„ » eg, with 



and 



«£ = s(«K + «K-.)±sa/(«£-«S-.) + ! W (25) 



45 = 5 (4!,! + i + <£) =f y (4!>« - + (26) 

As in the above, the signs should be chosen so that eft — > eft and eft — > eft for Q c — > 0. 
One of the resulting quasienergy curves is represented by crosses in Fig. E(a) • As seen from 
the figure, it reproduces the exact quasienergies well. The only significant differences occur 

(2) (2) 

when A c w A C) ± m with m > 2, where the approximate quasienergy manifolds e^ and e c ,n 
cross rather than anticross each other. The approximation can be improved by taking into 
account more off-diagonal components of the off-diagonal blocks of H' 2F and further iterating 
the process leading to Eqs. ( 12 5 p and The effective Rabi frequency for two quasienergy 
curves intersecting at A c ps A Cj ± m is Q c | J± m (E c /2cj T {)\. The gap of the corresponding avoided 
crossing is therefore a decreasing function of m when, as we have been assuming throughout 
this discussion, |E c /(2o; r f)| 1. 

Treating the crossings at A c w A C) ±i separately from that at A c = E c is justified if the 
former are sufficiently far apart from the latter, namely if |A Cj ±! — E c | ^> max (|f2 c J |, |O c Ji|). 
Given that Jo(E c /2a; r f) ^> Ji(E c /2a; r f) when E c /(2a; r f) 1, this condition can also be 
written 

f2 c Jo 



(27 > 

The iterative method outlined above fails when the inequality (|27|) is violated. The crossings 
between the manifolds may then occur at different values of A c than predicted by Eqs. f f23|) 
and (EU). 

Finally, we note that for weak fields the quasienergies obtained by the iterative method 
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are consistent with those predicted by perturbation theory, namely, for n = 36] 

0W2) 2 



e b,0 



e c,0 



E 



X c — A c — 2nu 



if 



(28a) 
(28b) 



3. The N -level approximation 

As seen above, the bare state \b) turns into a manifold of sideband states when the system 
is described in the Floquet formalism. The different states of this manifold are not directly 
coupled to each other by the rf field under our assumption that the dressing of state \b) 
by this field is negligible; however, they interact with each other indirectly, through their 
coupling with the harmonic components of state |c) by the control laser field. Neglecting 
the manifold structure of state \b) amounts to setting c& >n = for n ^ in Eq. fISaj) and to 
reducing the Hamiltonian matrix H'2F to the matrix 



/ 



2F 





ficJ-i/2 
^c</+i/2 



V 



fi c j_i/2 n c j /2 n c J+i/2 

-A c + 2w rf 
S c -A c 

S c - A c - 2w rf • • 

: ; ; ■.} 



(29) 



We refer to this approximation as the iV-level approximation, in view of the fact that the 
matrix H^p is effectively the Hamiltonian of a system consisting of a single state \b', 0) = \b) 
interacting with a manifold of N independent states |c', n)), N — > oo — see Fig. \U(b). As 
mentioned in the Introduction, this approximation has often been made in studies of the 
interaction of low-lying states with rf-dressed Rydberg states. 

The corresponding quasienergies are the eigenvalues of the matrix H^p. An example of 
the resulting spectrum is shown in Fig. &[b). Proceeding as in Sec. Ill B 21 one finds that for 



tt c = the corresponding quasienergies are equal either to e 



~(o) 



Or tO €c°n 



(0) c 

e c .n for some 



value of n, whereas for Q c ^ and |E c /(2a; r f)| 1, they are approximately equal either to 



-c,0 



e^d, or to 6c,n = £ { c,h(n ^ 0). In the latter case, the c-manifold intersects 



.(0) 



the fe-quasienergy curve at values of A c for which e^ 1 ^ = ej,n for some value of n. However, 



13 



these intersections normally occur at different detunings than found in Sec. Ill B 21 since in 
general ei,n 7^ e^n when n^O. Instead, the crossings occur at A c = A Ci ± m where 

A Cj±m = S c ± 2mu rf =f Q 2 c Jq I (8mu rf ). (30) 

Comparing this result with Eq. (1241) , we see that within the iV-level approximation the shift 
of the crossings from their zero-f2 c positions is about half the correct value [37 ] . 



C. Dressed state dynamics and relaxation 

Being entirely based on Eq. ([6]), the analysis developed in Sec. Ill Bl neglects both the 
frequency width of the control laser and the energy width of the states \b) and |c). The 
corresponding time evolution of the system is purely oscillatory, with population transferred 
to-and-fro between \b) and |c) without any damping. The evolution would be a pure Rabi 
flopping of characteristic angular frequency (A 2 + f^ 2 ) 1 / 2 in the absence of the rf field. The 
oscillation is multimode in the presence of the latter, with Fourier components separated in 
angular frequency by integer multiples of 2cu r f. Eqs. fl2Tj) and (122]) approximately give the 
corresponding characteristic frequencies as [(E c — A c ) 2 + f2 2 Jq} 1 ^ 2 /2 + 2ptu r f, p = 0, ±1, ±2, . . . 

Relaxation modifies this picture. Taking into account the natural width of state |c), T c , 
and the linewidth of the control laser, 7 C , and assuming that |c) can de-excite only to \b), 
the time evolution of the system is described by the following optical Bloch equations: 

Pbb = r c p cc + -^-{pbc~ Pcb), (31a) 

Pcc = -r c p cc + -^-(pd,- Pbc), (31b) 

Pbc = ~ {r c /2 + 7 C + i [A c - 2S C sin 2 u r{ t] } p bc 

+ ~Y(pbb- Pec)- (31c) 

(Recall that at the moment we study the system formed by states \b) and |c) in isolation 
from state \a). We therefore assume, for the time being, that \b) does not de-excite to a 
lower state. The finite lifetime of state \b) will be taken into account in the analysis of the 
full 3- level ladder system, in Sec. IIIII ) Decoherence dampens the dynamical evolution of the 
system until the oscillation is purely driven by the rf field. In this steady state, which is 
reached for t ^> 1/T C , the elements of the density matrix oscillate at integer multiples of the 
fundamental angular frequency 2w r f about a constant cycle-average. 
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In the iV-level approximation, we treat each side band of state |c) as an independent 
bound state and describe the evolution of the system by optical Bloch equations based on 
the Hamiltonian matrix H' 2F , namely 

Pbb = r c ^Pfcfc + ^(JkPbk - JkPkb), (32a) 
k k 

Pbk = - [r c /2 + 7 C - i(E c - A c - 2kua)] pbk 




(32b) 



p jk = - [r c - 2i(j - k)u Ti ]p jk + -^(JkPjb ~ JjPbk), (32c) 

with the indexes j and k running over all the states forming the c-manifold. In this approach, 
the population in state \c), p cc , is defined as 1 — p^, or equivalently, as p cc = J2 k p k k- In 
contrast to Eqs. ()31aH31cj) . which predict oscillatory populations and coherences, Eqs. f)32af - 
I32c[) lead to a constant density matrix in the long time limit. 

We assume that the atom is initially in state \b) and calculate the population transferred 
to state |c), p cc (t), and its temporal average in the steady state, (p cc (oo)). We define the 
latter as 

1 f t+T 

(Pcc(oo)) = - lim / Pcc (t')dt', (33) 

with T = 2ir/ (2o; r f). We thus set p^ to 1 and all the other elements of the density matrix to 
zero at time t = and numerically solve either Eqs. fl31aH31c| or Eqs. f)32aH32c| from t = 

to t > i/r c . 

How the population transferred to state |c) varies with the detuning A c depends both 
on the value of T c and on whether the inequality ( 1271) is fulfilled or not. Avoided crossings 
between quasienergy curves can be treated in isolation when this inequality is fulfilled, and 
adjacent resonances do not overlap significantly due to the natural width of the upper level 
when T c 2u T f. The variation of (p cc (oo)) in these conditions is examplified by Fig. EJ^c). 
The results represented respectively by a dark pink curve and by a light purple curve were 
obtained by solving Eqs. fl31aH31cj) and Eqs. f l32aH32cj) . 

As seen from the figure, the avoided crossings in the dressed state spectrum are associated 
with strong enhancements of population transfer. The origin of these enhancements is easily 
understood given that T c and 7 C are both sufficiently small that the evolution of the system 
is dominated by its dressed states dynamics rather than by relaxation. For the parameters 
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considered, the dressed atom can be described as a simple two-state system in the vicinity 
of each avoided crossing. Sufficiently far from any crossing, the dressed states are close to 
either one of the \b',n)) states or one of the \c',n)) states, with little admixture between 
the two. At a crossing, however, the interacting dressed states are approximately equal 
superpositions of \b) and |c). Population transfer is therefore enhanced at each crossing, 
typically over the range of values of A c for which the relevant b- and c-dressed states are 
significantly admixed with each other. Treating the dressed atom as a pure two-state system 



in the vicinity of each crossing leads to a simp 



to the upper state, i.e., for A c pa E c + 2nu T { 



e expression for the population transferred 



Q 2 c Jl/2 

(Pcc(oo)) ra r c (r c /2 + 7c ) [i + (A c - i c - 2^)2/^] + opf (34) 

Thus (p cc (oo)) exhibits power-broadened sidebands whose height and width are modulated 
by the square of the Bessel function factors J n . When \n\ is larger than the argument of 
these Bessel functions, these enhancements decrease in amplitude and their width tends to 
the natural width of the upper level, T c , for increasing sideband orders. For the parameters 
of Fig. [2(c), this analysis predicts that the n = 0, n = ±1 and n = ±2 enhancements have 
an heigth of 0.50, 0.46 and 0.08, respectively, and a full width at half maximum of 13.3, 
3.6 and 1.1 MHz. As can be seen from the figure (dark purple curve), these values are in 
agreement with the numerical results, although the n = sideband somewhat overlaps the 
others. (The inequality (]2"T]) is fulfilled but marginally here, since \Q C Jo\/(2u r {) « 0.5.) 

The picture is the same in the iV-level approximation for the parameters considered. The 
discrepancy in the position of the enhancements mirrors the differences in the position of 
the crossings in the quasienergy maps. It can be largely eliminated by a simple fc-dependent 
shift in the S c — A c — 2ku r { factor multiplying p bk in Eq. (I32b[) j^. 

The comparison shows that for the parameters considered in Fig. [2} the upper state ef- 
fectively turns into a manifold of states well described within the A/"- level approximation. 
As mentioned in Sec. Ill B 2\ increasing Q c eventually results in a break down of this approx- 
imation as the states forming the Floquet manifold will start to interact with each other. 
For sufficiently strong control fields, the resonance sidebands overlap each other and merge 
with the zeroth order resonance. 

How large \Q C Jq\ is relative to 2cu r f also impacts on the temporal evolution of the system. 
When |f2 c Jo| <C 2u; r f, the atom interacts with the laser field as if state |c) were a manifold 
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FIG. 3: (color online) The population in the upper state at sin 2 u; r fi = (first row), 1/2 (second 
row) and 1 (third row), and its time average (bottom row), in the steady state regime (t 3> 1/T C ). 
£ c /(2vr) = 20 MHz and 7 C = in all cases, (a): Q c /2tt = 20 MHz, r c /2vr = 1 MHz and w rf /27r = 5 
MHz. (b): n c /2ir = 0.1 MHz, r c /2vr = 1 MHz and w rf /27r = 5 MHz. (c): O c /2vr = 0.1 MHz, 
r c /2vr = 10 MHz and w rf /27r = 5 MHz. (d): O c /2vr = 0.1 MHz, r c /2vr = 10 MHz and w rf /2vr = 0.01 
MHz. S c /(27r) = 20 MHz and 7 C = in all cases. p cc is multiplied by a factor of 0.025 in column 
(b). 



xO.025 




of stationary states. In the opposite limit, the rf field changes slowly on the time scale on 
which the states \b) and |c) interact with each other. This interaction then occurs as if the 
rf field was static at any instant. In this case, the Stark shift due to the rf field is expected 
to make the population in the state |c) oscillate periodically. 

Columns (a) and (b) of Fig. [3] compare instantaneous and cycle-average values of p cc in 
the steady state regime for two different values of Q c . In column (a), |f2 c J |/(2o; r f) ~ 0.4. 
This value is insufficiently small for the iV-state picture to hold well, which leads to a clear 
oscillation in the population in the state |c). In column (b), |f2 c J |/(2o; r f) ~ 0.002, which 
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fulfills the condition ( 135]) : p cc is almost constant in time. The resonance sidebands are also 
much better marked than in column (a) because of the much smaller power broadening. 

The population in the state |c) also depends on the time scale of the relaxation mecha- 
nisms. Well defined enhancement sidebands are not expected unless spontaneous decay is 



slow compared to the oscillation of the rf field, i.e., unless 13 1 



T c < 2co r{ . (35) 

This inequality also expresses the requirement that adjacent sidebands, which are separated 
roughly by 2oj r f, do not overlap due to the the natural linewidth of the resonance, T c . The 
role of the lifetime of the state |c) is illustrated by column (c) of Fig. [31 which was calculated 
for the same parameters as column (b) apart for T c : while r c <C 2u r f in column (b), r c = 2u; r f 
in column (c). In the latter case the sidebands partly overlap each other and the population 
in the state |c) oscillates markedly. 

Finally, column (d) of Fig. [3] shows the time evolution of p cc for the same parameters as in 
column (c) apart that 2u r { is much smaller than both r c and A c . The N- level picture is now 
invalid and there is no sign of sidebands. Instead, the profile of p cc is a single peak of width 
r c following the instantaneous position of state |c) as this state is periodically Stark-shifted 



by the rf field 13| 



III. RF-DRESSED THREE-LEVEL SYSTEMS 



A. General formulation 



We now take into account the coupling of state \b) to state \a) by the probe laser field and 
allow state \b) to de-excite to state \a), which is assumed to be stable (Fig. H]). We proceed 
as in Sec. Ill Al and, ignoring relaxation at first, write the state vector of the atom as 



Mt))=c a (t)\a)+c b (t)\b) + c c (t)\c), 



(36) 



with the amplitudes c a (t), Cb(t), c c {t) satisfying the equation 



d_ 

dt 



V c c/ 



-zfi p /2 
-itt p /2 iA p -iQ c /2 
, —ifl c /2 «[Ar — 2S C sin 2 w r ft]y 



Cb 

V c c/ 



(37) 
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FIG. 4: (a) Schematic of the rf-dressed three-level Rydberg system. The probe laser (angular 
frequency w p ) couples the states \a) and while the control laser (angular frequency u) c ) couples 
the states \b) and |c). The energy of the latter is modulated by the rf field (angular frequency w r f). 
The detunings of these two laser fields from the respective transition frequencies are A p and A c . 
The state |c) decays to the state \b), and the state \b) to the state \a). (b) The manifold structure 
of the states \b) and |c) induced by the rf field. The energy separation between adjacent sidebands 
is 2u> r f. The rf field also shifts the energy of state |c) by an amount S c ; the effective detuning of 
the control field is thus A c — E c . 



Following the same analysis as in Sec. Ill we formulate the problem in terms of the Floquet 
Hamiltonian 



where 



and 



( a n p o ^ 



H 



3F 



fit B fi c 
fit C 



A = diag(. . . , 2cj t {, 0, — 2cj t {, . 



fi r 



^~^p ^~^p ^"^p 



diag(..., T , T , T , 




H 



2F 



1 
1 



(38) 



(39a) 
(39b) 

(40) 
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Generalizing the discussion of last section, the corresponding Floquet states are described 
by time-independent vectors \^k)) which can be written in terms of the vectors \a, n)) = 
(\a)5 in ,i = 0, =Fl, =f2, . . .) T and of the vectors \b,n)) and \c,n)) defined by Eqs. (TTHT) and 
(ED: 

oo 

I**)) = ^2 (c a -n\a,n)) + c b _ n \b,n)) + c c _ n \c,n))). (41) 

n=— oo 

As before, the C-block of H3F can be diagonalized by a change of basis transforming the 
matrix H 2 f into the matrix H' 2F defined by Eq. ( flT|) . The change transforms the matrix H 3 p 
into the matrix H 3F differing from H3F by the replacement of H2F by H^p. It does not affect 
the spectrum of quasienergy. Making the iV-level approximation further reduces H 3F to 



ft 



3F 



/ fi p /2 

Q p /2 -A p 

; 

fi c J_i/2 

rt c j /2 

n c j+i/2 



000 

Vl c J_\/2 Vt c J§/2 Vt c Jj r \/2 



A R + 2u r{ 
S c - A R 

S r 






Ar — 2u; r f • 



(42) 
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The optical Bloch equations for the 3-state system with rf-modulation can be written as 

Paa=TbPbb + -^-{.Pab- Pba) , (43a) 
Pbb = — FbPbb + r c p cc -^-(Pab — Pba) 

+ ^(p bc -Pcb), (43b) 

Pcc = -^cPcc T^(Pbc-Pcb), (43c) 

Pab = ~ (^A p + — + 7 p )p a6 

- ~Y{pbb ~ Paa) + -^Pac, (43d) 

p fec = - (z'A c - 2z'£ c sin 2 w rf t H — - + 7 P + 7c)p& c 

- ^ (Pec -Pbb)~ ^Pac, (43e) 

r 

Pac = - (^A p + iA c - 2iS c sin 2 u v{ t + y + lc)p ac 

+ ^P ao - ^p 6c - (43f) 

As in Sec. HI1 T c and 7 C denote the natural width of state |c) and the linewidth of the control 
laser. 1^ and 7 P , which were not introduced in Sec. [Til are t ne corresponding quantities for 
state \b) and the probe laser (Fig. HJ). Correspondingly, the optical Bloch equations in the 
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iV-level approximation now read 



Paa — 1 6Pb6 H 7^{Pab ~ Pba) j 



P&b 



FbPbb 



(Pab - Pba) 



+ r c ^ p fcfe + — ^ ^( JfcP^ - JfcPfc ), 



— + 7p + 2Ap p ab + —{paa ~ Pbb) 



(44a) 



(44b) 



(44c) 



Pak — — 



+ 7 P + 7c + ^(A p + A c - £ c + 2ku- 



rfJ 



Pafc 



Pbk 



2 "Pfcfc 4 



"Pab, 



(44d) 



+ 



Pbb 



+ 7c + «(A c 



E c + 2/cw r f) 



Pbk 



'Pak, 



Pjk 



[r c - 2z(j - A;)w r f]p jfc + -^-{JkPjb ~ JjPbk), 



(44e) 
(44f) 



with the indexes j and k running over all the components of the Floquet manifold of state |c). 
These equations are to the Hamiltonian H' 3F what Eqs. f l43aU43ff are to the Hamiltonians 
H 3 f and H' 3F . As in the 2-state case, the solutions of Eqs. fl44aH44f|) tend to constants in the 
long time limit, whereas those of Eqs. ( l43aH43ff tend to a steady state in which the elements 
of the density matrix oscillate at a fundamental angular frequency 2u T f. In either case, the 
absorption coefficient for the probe beam is given by the equation 



a 



-k p Im 



27rN|d ba • e F 



(P6a(00)) £ 



(45) 



e htt p 

where k p is the wave number of the probe field, N is the number density of the atoms forming 
the medium, db a is the matrix element of the dipole operator between the states \a) and 
| b), e p is the polarization vector of the probe field, e is the permittivity of free space, and 
(pb a (oo)) av denotes the average value of p ba in the stationary regime. For a cold atomic 
ensemble, (p ba (oo)) av = (p ba (oo)), with 

1 f t+T 

(p 6a (oo)) = - lim / Pba {f)'dt'. (46) 
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For a thermal ensemble where the Doppler effect must be taken into account, (pf, a (oo)) av is 
the average of (pf, a (oo)) over the Maxwellian distribution of velocity of the atoms. 

B. Absorption sidebands and induced sidebands 

1 . Weak control fields 

A representative section of the quasienergy spectrum of the full Floquet Hamiltonian 
is shown in Fig. 15(a) for a case where the control laser is sufficiently weak that the in- 
equality (I27p is fulfilled. The horizontal line interrupted by avoided crossings indicates the 
quasienergy of the zeroth order Floquet sideband of state \a). The oblique lines arise either 
from the sidebands of the state \b) or from those of the state |c), as indicated in the caption. 
The b- and c-quasienergy manifolds run parallel to each other in Fig. [5(a), although they 
intersect in Fig. [2(a) , because the parameter varied in Fig. [5] is the detuning of the probe 
field while that varied in Fig. [2] is the detuning of the control field. The latter is actually 
set to zero in these calculations, which means, in view of the Stark shift of state |c), that 
the control laser is off resonance. The quasienergy spectrum obtained in the A-level ap- 
proximation is shown in Fig. [5(b) , for comparison. The absorption profile for a cold atom 
ensemble is given in Fig. [5(c) . (The quantity represented in Fig. [5(c) is (p;, a (oo)), which 
is proportional to the absorption coefficient if the atomic motion can be neglected. The 
maxima of (pb a {oo)) correspond to maxima of absorption.) 

The dressed states crossing the horizontal lines all contain the bare state \b) — i.e., at 
least one of the c^ n coefficients is non-zero in Eq. fHTj) . Therefore each of these crossings 
corresponds to a resonant coupling between this state and the bare state \a) and may give 
rise to an enhancement of absorption of the probe field coupling \a) to \b). Fig. [5(c) shows 
a large enhancement of absorption at the crossing involving the zeroth order sideband of 
state | b) (at A p = 0), and smaller enhancements at the other crossings. The enhancement 
at the crossings with the +1 sideband of |c) (at 85 MHz) is hardly visible, though, and the 
enhancements associated with the crossings of the +1 and +2 sidebands of \b) (at 50 and 
100 MHz, respectively) are too small to be seen. These enhancements are weaker compared 
to those involving the zeroth-order sidebands both because of the scaling of fl c by the J n 
factors ( | J n | is much smaller than 1 when n^O for the parameters of Fig. [5]) and because 
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A p /2tt (MHz) 

FIG. 5: (color online) The quasienergy spectrum, (a) and (b), and the absorption profile of the 
probe field, (c), in an rf-dressed 3-level system. O p /27r = 1 MHz, Q c /2ir = 8 MHz, T, c /2tt = 35 
MHz, A c /2vr = MHz, w rf /2vr = 25 MHz, r fe /2yr = 6 MHz, r c /2vr = 0.01 MHz and j p /2it = 
7c/27r = 0.1 MHz. (a) and thick dark blue curve in (c): exact results, (b) and thin light pink curve 
in (c): results in the iV-level approximation. In (a), the crosses identify e^o and the circles e C) o, e c ,l 
and e c -\. The other oblique lines correspond to e^^, n ^ 0. The horizontal line corresponds to e a> o- 
One may note small differences in the position of the crossings between (a) and (b). Because of 
the significant Stark shift of the upper state, the control laser is off resonance although A c /27r = 0, 
which explains why the absorption profile of the probe field does not exhibit EIT dips. 

of the larger energy difference between the respective sidebands (a larger energy difference 
means a smaller admixture of \b) into the sidebands of \c)). However, as will be seen shortly, 
significant enhancements may occur at the crossings with the n ^ sidebands of b for 
stronger control fields. 

Comparing Fig. [5^a) to Fig. [5^b) and the thick blue line to the thin pink line in Fig. 
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EJc), we see that the N- level approximation holds well in the present case, as could have 
been expected since the condition flM|) is fulfilled. Because the side bands of the 6-state 
are neglected in the iV-level approximation, the 6-quasienergy manifold of Fig. E{a) reduces 
to a single quasienergy curve in Fig. |^b) . However, the c-manifold is practically the same 
in the two spectra. The absorption profile is also very well reproduced in the iV-level 
approximation for the parameters considered: the two sets of results are almost identical 
apart for very small differences in the position of the enhancements (noticeable only for the 
enhancement at A p = — 15 MHz on the scale of the figure) and very small differences in 
their amplitude. 

2. Strong control fields 

We now increase Q c /2n to 30 and 60 MHz, keeping the other parameters the same as 
in Fig. [51 The resulting absorption profiles are shown in Figs. E](a) and[6](b), respectively. 
As expected from the previous discussion, the iV-level approximation deteriorates as the 
strength of the control field increases. It is still in rough agreement with the exact result 
in the case of Fig. [6]^a), for which Q c is smaller, but not much smaller, than 2w r f, although 
there are differences in the position and the amplitudes of the enhancements. However, 
when the condition (13"1"|) is more strongly violated, the iV-level approximation breaks down 
completely. There is little agreement with the exact result in Fig. [f^b). 

An interesting feature of Fig. [6](a) is the presence of a small enhancement of absorption 
at A p /27T « 48 MHz. This enhancement is not present in the iV-level approximation. It 
coincides with the crossing of \a) with the +1 sideband of \b) in the quasienergy spectrum, 
at 2hu r { above the main absorption peak (which is concomitant with the crossing of the 
0-th order sideband of \b).) This feature is more prominent in Fig. MJo), at A p /2tt 30 
MHz; the +2 sideband of \b) is also (barely) visible, at about 80 MHz at this higher power 
of the control laser. We refer to such sidebands as "induced sidebands" , as they arise from 
sidebands of \b) induced by the coupling of this state with the sidebands of |c) by the control 
field. (Recall that the rf field is assumed to be too weak to dress state \b) directly: the b- 
manifold originates from the coupling of \b) to the c-manifold by the control field, not from 
the coupling of \b) to itself by the rf field.) 

That the sidebands of \a) play no role in this absorption profile is shown by the good 
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agreement between the exact results (thick blue curves) and the results obtained when these 
sidebands are neglected altogether (green dashed curves). 

Before closing, we briefly comment on whether the conditions Eq. ( 1271) and ( 1351) are 
sufficient to garantee the validity of the iV-level approximation for a three-state system. 
Recall that these conditions mean that the oscillations of the rf field are faster than the 
time scales over which the system evolves under the effect of the control laser field and the 
spontaneous decay of state |c). However, the probe field introduces an additional time scale 
in the problem: if this field is excessively strong, it will disturb the dynamics of the b-c 
system too rapidly for the state |c) to show a manifold structure. In this case, the Floquet 
sidebands of state \a) induced by its coupling to state |c) can no longer be ignored. Given 
that the effective Rabi frequency for the coupling of \a) to |c) is approximately given by 
^p^c/2S c 
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39l | . we expect that the condition 



l^l«2c rf (47) 

should be added to the conditions Eq. (1271) and fl35l) for the iV-level approximation to be 
valid. Should the inequality ( H71) not be fulfilled, then it is likely that the Floquet sidebands 
of | a) need to be taken into account, as well as those of \b) and |c). 



C. EIT in rf-dressed atomic ensembles 



1. Cold atomic ensembles 



We still assume that the Doppler effect is negligible, as in the last section. However, we 
now set A p to and vary A c rather than set A c to and vary A p . If the rf field was absent, 
the control laser field would couple the states \b) and |c) resonantly at A c = 0, which would 
be accompanied by a dip in Im[(p a b(oo))]) . (This is the well known EIT absorption window 
at A p = In the presence of the rf field, this dip is displaced to A c = S c , due to the 
ac Stark shift of |c), and splits into multiple sidebands. 

EIT with rf modulation is illustrated by Fig. [7J These results were calculated for the 
same values of fl p , S c , cj t {, Tf,, T c , 7 P and 7 C as in Figs. [3]and[5l The probe Rabi frequency, 
Q p /2n, increases from 8 MHz in Fig. 0(a) to 30 MHz in Fig. [TJb) and to 60 MHz in Fig. 
0(c). As seen from the figure, the rf field creates transmission sidebands on either side of 
the main EIT feature. Since A p is set to zero and the frequency of the control field is varied, 
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FIG. 6: (color online) The probe absorption in the same rf-dressed 3-level system as in Fig. [5] but 
for stronger control fields for which the N- level approximation is inaccurate. Q, c /2tt = 30 MHz in 
(a) and 60 MHz in (b). Thick blue solid curves: results obtained by solving the full optical Bloch 
equations, Eqs. (|43aH43fl) . Thin pink solid curves: results obtained in the iV-level approximation, 
Eqs. (|44aH44f|) . Dashed green curves: results obtained by adding the sidebands of state \b), but 
not those of state \a), to the N- level approximation. The enhancements associated with crossings 
between the lower state and the intermediate state manifold are identified by the letter b between 
the lower state and the Rydberg state manifold by the letter c. 



these sidebands occur at the values of A c at which the c-quasienergy manifold intersects the 
zeroth-order quasienergy sideband of the state \b). In the weak control case of Fig. [7(a), 
they manifest as narrow dips regularly spaced by 2u; r f 40(. Apart for minor differences in 
the depth, the iV-state approximation reproduces these features very well (the thick blue 
curve to the thin pink curve). The absorption profile is thus well explained by the model in 
which the upper state simply splits into a comb of Floquet states, each one interacting with 
the intermediate state independently from the others. However, this picture changes when 
Q c approaches or exceeds 2u T { — see, respectively, Fig. [7(b) and (c): the EIT dips broaden 
and shift towards the zeroth-order sideband when Q c increases, and the absorption profile 
first acquires a peak-and-trough structure and then changes into a series of absorption peaks 
(rather than absorption dips) superimposed on a slowly varying background. These peaks 
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A c /2tt (MHz) 

FIG. 7: (color online) Probe absorption vs. control laser detuning in an rf-dressed 3-level system, 
for three different control Rabi frequencies. The parameters of the system are the same as in Figs. 
[5] and El but here the probe field is on-resonance (A p = 0) and A c varies. Q c /2tt is 8 MHz in (a), 
30 MHz in (b) and 60 MHz in (c). Thick blue solid curves: results obtained by solving the full 
optical Bloch equations, Eqs. (|43aH43f|) . Thin pink solid curves: results obtained in the iV-level 
approximation, Eqs. (|44al I44fj) . Dotted green curves: absorption profile without rf field (the energy 
of the Rydberg state was shifted by S c to facilitate the comparison with the other results). As 
Q c increases, the EIT transparency window in (a) expands into the sidebands region and the EIT 
absorption dips change into absorption peaks. The width of the features is modulated by the 
square of the Bessel function factors J n . 

occur at the detunings at which state \a) is resonant with the dressed state formed by the 
states | b) and |c) coupled by the control field. The same changes are observed in the N- level 
approximation, at least for the parameters considered, but the position of the absorption 
features is more and more incorrect. 

The green dotted curves represent the absorption profile calculated without the sidebands 
of state |c) (but taking its Stark shift into account). Comparing these results to the thick 
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FIG. 8: (color online) The probe absorption in the same rf-dressed 3-level system as in Fig. [7] 
but here with Doppler averaging and plotted against probe detuning. The probe and control laser 
beams are assumed to be colinear and counterpropagating. S7 c /27r is 8 MHz in (a) and 30 MHz in 
(b). A c = in both graphs. Thick green curves: results obtained by solving the full optical Bloch 
equations, Eqs. (|43aH43f|) . Thin black curves: results obtained in the iV-level approximation, Eqs. 
(|33aH44i|) . The higher power of the control field in (b) results in a broadening and a shift of the 
EIT dips and to inaccuracies in the iV-level approximation. 

blue curve in Fig. [7(c) shows that the suppression of EIT in narrow ranges of frequencies at 
large control powers arises from the interplay between these sidebands. 

2. Thermal atomic ensembles 

Let us assume that the two laser beams are colinear and counterpropagating. Compared 
to an atom at rest, the probe and control angular frequencies for an atom moving with 
a velocity component v in the direction of the control laser beam are shifted respectively 
by k p v and —k c v, where k p and k c are the wave numbers in the laboratory frame. In the 
reference frame of the atom, the probe detuning is therefore upshifted to A p + k p v and 
the control detuning downshifted to A c — k c v, where A p and A c are the detunings in the 



29 



laboratory frame. The absorption of the probe beam at A p and A c will thus be determined 
by the velocity average of the steady-state coherence (p a &(oo)) calculated at the Doppler- 
shifted detunings. We denote this average by (pjj(A p , A c , oo)) av . For atoms of mass m at a 
temperature T, 

/oo 
/(f)(py(A p + k p v, A c - k c v, oo)) dv (48) 
-oo 

with 

where k-& is the Boltzmann constant. 

Fig. [8] shows the probe absorption calculated using the A^- level approximation (thin black 
curves) or calculated directly from Eqs. f l43aH43f]) (thick green curves), as a function of the 
probe detuning, for two different values of Q c . The other parameters are the same as those 
used in Fig. |5j In the Doppler average, we assume a cloud of 85 Rb atoms at a temperature 
of 40 °C, and probe and control wavelengths of 780 and 480 nm, respectively. 

We see that the TV-level approximation is accurate in the case of figure (a), for which the 
control power is relatively low. For small values of Q c , each absorption minimum corresponds 
to EIT in the velocity class which simultaneously satisfies the two conditions A p + k p v = 
and A c — k c v — £ c + 2nu v { = 0. For a fixed value of A c , the corresponding absorption dips 
occur at the probe detunings 

A p = ^(S c - A c -2nu li ). (50) 

Adjacent EIT dips are thus separated in probe detuning by 2u) T fk p /k c in these conditions. 

Increasing Q c both widens the EIP dips through power broadening and shifts their posi- 
tions. The latter effect is due to the shift in the sidebands of |c) arising from their coupling 
with state \b). When this shift is non- negligible, the second of the above resonance condi- 
tions must be replaced by A c — k c v — S c + 2nu r { — 5 n = 0, where 5 n is an f2 c -dependent shift, 
and the absorption dips occur at the probe detunings 

A p = ^(S c - A c -2nuj rl + 5 n ). (51) 

As seen from the figure, the A^-level approximation becomes inaccurate when Q c / (2tt) is 
increased from 8 to 30 MHz: at the higher control power there is a clear difference both 
in the position and in the depth of the EIT dips between the absorption profiles calculated 
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with and without making this approximation. (Correcting the TV-level approximation for 
the sidebands of the intermediate state, still neglecting the sidebands of the lower states, 
restores the agreement with the exact results.) 



IV. SUMMARY 

To conclude, we have investigated the absorption spectrum of two- and three-level model 
atomic systems with ac-modulated quadratic Stark shift, both in the weak coupling regime 
and beyond this regime. Although we specifically consider the case of a Rydberg system 
modulated by a radio-frequency field, our results are generic. 

The different regimes investigated can be defined in terms of the modulation amplitude 
of the upper state energy, which is characterized by the quantity S c defined by Eq. ([7j), and 
in terms of the angular frequency of the rf field, w r f, of the Rabi frequency of the control 
field, Q c , and of the natural width of the upper state, T c . When the conditions |E C | 2a> r f, 
|O c J (S c /2a; r f)| <C 2w rf and T c <C 2u r { are all met, the quasienergy curve crossings in the 
Floquet spectrum (and the associated structures in the absorption spectrum) can be treated 
in isolation of each other. The system then behaves as if the most polarizable state is 
effectively a manifold of sideband states and the approximation we refer to as the TV-level 
approximation is accurate. This approximation amounts to neglecting the manifold structure 
of the less polarizable states and model the ac modulated system as a many-level system in 
which these states interact with a manifold of independent sideband states spawned by the 
most polarizable state. 

Increasing the strength of the control field to |f2 c Jo(E c /2o; r f)| ps 2u r { shifts the position of 
the absorption sidebands and changes their amplitude, in agreement with the perturbative 



analysis developed in a different context in Refs. [25J and [27| and with the iterative approach 
outlined in Sec. Ill B 2[ The A- level approximation becomes inaccurate as \Q C \ increases, in 
particular in predicting sidebands shifts half too small. However, this approximation can be 
brought into agreement with the correct positions of the sidebands by a simple change in 
the optical Bloch equations (or the corresponding Hamiltonian) . 

In the strong coupling regime, where |f2 c J (E c /2o; I .f)| ^> 2u r f, the A-level approximation 
fails unless it is corrected by including the intermediate-state Floquet manifold in the model. 
For sufficiently large values of |fi c |, the probe absorption spectrum in a 3- level ladder system 
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develops sidebands induced by the coupling of the intermediate state to the polarizable upper 
state, in addition to the sidebands manifesting in the weak coupling regime. For certain 
combinations of the parameters, this spectrum changes from one exhibiting sharp EIT dips 
to one exhibiting sharp absorption peaks as \Q C \ increases from |f2 c J (E c /2o; r f)| <C 2w r f to 
|fi c J (S c /2w rf )| > 2w rf . 
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